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Abstract 



Using a quantum formulation of the master equation we study a kinetic 
Ising model with competing stochastic processes: the Glauber dynamics with 
probability p and the Kawasaki dynamics with probability 1 — p. Introducing 
explicitely the coupling to a heat bath and the mutual static interaction of the 
spins the model can be traced back exactly to a Ginzburg Landau functional 
when the interaction is of long range order. The dependence of the correla- 
tion length on the temperature and on the probability p is calculated. In case 
that the spins are subject to flip processes the correlation length disappears 
for each finite temperature. In the exchange dominated case the system is 
strongly correlated for each temperature. 
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I. INTRODUCTION 



The time evolution of nonequilibrium Isine. spin systems is of interest in particular when 
the spins are subject to different dynamicslH. Initiated by Glauberl, for a review seel, a 
stochastic model had been considered in which each spin can flip with a certain transition 
rate constructed in such a way that detailed balance is fulfilled. The system has to be in 
contact with a heat bath not specified in detail. The bath yields the energy for such a 
stochastic flip-process where every spin can flip in a first view independently on the orienta- 
tion of neighbouring spins. In general it is assumed that the dynamics of the whole system 
is governed by a stochastic master equation. Recently, a very powerful method in analysing 
master equations was proposediH3. It based on a quantum formulation of the underlying 
master equation written in terms of creation and annhiliation operators. Originally, the 
method had been formulated for Bose operators3H§. More appropriate, in particular for 
problems with an exclusive dynamics, is the application of quantum Pauli operators which 
allows exact solutions for a whole class of one dimensional systems, for a review see @. We 
have proposed an alternative way to include the coupling to a heat bathS explicitely. The 
method has been applied to study a spin-facilitated mode!0 where the dynamics is re- 
stricted by specific constraints. It can be also used to find the above mentioned transition 
probability^ which is not unique. It should be emphasized that the method is determined 
by the algebraic properties of the operators where an extension to a p-state modelEl and to 
a Q-statistics is straightforward^. 

Using such a Fock space formulation of the master equation we analyse a model where the 
system experiments with propability p single spin flips (Glauber dynamics!) and with propa- 
bility 1 — p spin exchange processes (Kawasaki dynamics!). The coupling to a heat bath 
with fixed temperatureT is included and the static interaction between the spins is assumed 
to be of infinite ranged. Whereas the Glauber kinetics is always related to a change of the 
order parameter (nonconserved order parameter) , the Kawasaki dynamics simulates the flux 
of energy into the system (conserved order parameter). In a continuous approach the model 
consists of a combination of model A and model B within the classification proposed by 
Hohenberg and Halperin@. 

There is a great effort in analysing such nonequilibrium spin systems with competing dy- 
namics. Gonzalez-Miranda et al£3 had studied a kinetic Ising model where Glauber and 
Kawasaki dynamics drive the system simultaneously. Using Monte Carlo simulations in 
two dimensions they found the phase diagram in the plane temperature versus the rate 
of Kawasaki process, here 1 — p. It reveals a line of continuous transitions between the 
ferro-and the paramagnetic phases. Assuming that the temperature depends on p there 
appears for low temperatures a nonequilibrium tricritical point which had been confirmed 
by DickmanEI employing the dynamic pair approximation. In the present paper, the rate p 
and the temperature T are assumed to be independent variables. 

Another interesting feature observed is the occurence of self organization phenomenal. In- 
creasing the energy flux the system changes continuously from the ferromagnetic to 
magnetic state. A further increase of flux drives the system into an antiferromagnetic phase 
which has been confirmed numerically at infinite temperature! and for finite temperatures!. 
Based on a Monte Carlo simulation Grandi et a!0 was sucessful in calculating the critical 
exponents when the Kawasaki dynamics dominats the behaviour of the system. An analyt- 
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ical approach has recently given by Ma et alii 

Here, we are interested in a quantum formulation of the problem using a Fock space repre- 
sentation of the master equation. The model can be attributed to the standard Ginzburg 
Landau functional exactly when the static interaction is realized by an infinite range model. 
Calculating the correlation function we find the correlation length depending on the inde- 
pendent variables temperature T and the rate p which is the probabilty that the Glauber 
dynamics is realized. In the model there is no indication for a tri critical point. 



II. FOCK SPACE REPRESENTATION 

The dynamics is govered by the master equation written in the symbolic form 

d t P(n,t) = L'P(n,t) (1) 

Here P is the probability that a certain configuration in terms of a lattice gas representation, 
n = {ni,n 2 . . .tin), is realized at time t. The evolution operator V specified below is given 
by competing flip and exchange processes. The variables are related to the occupation 
number operators with eigenvalues and 1. Hence, the problem is to formulate the dynamics 
in such a way that this constraint is taken into accountEiHll. The situation in mind can be 
analyzed in a seemingly compact form using a Fock space representation of the master 
equationilE introduced by eq. ([!]). Followingui00 the probability distribution P(n,t) is 
related to a state vector | F(t)) in a Fock-space according to P(n,t) = (n \ F(t)) where the 
basic-vectors | ft) are composed of second quantized operators. The master equation (HD can 
be transformed into an equivalent equation in a Fock-space 

d t | F(t)) = L | F(f)> (2) 

where the operator V in eq.(0) is mapped onto the operator L. Up to now the procedure is 
independent on the used operators. Originally, the method had been applied for the Bose 
caseQti Recently, an extension to restricted occupation numbers (two discrete orientations) 
was proposed!^!" based upon a Pauli-operator representation. These operators commute 
at different sites and anti-commute at the same lattice site. A further extension to an p-fold 
occupation number is possible^. 

The relation between the quantum-like formalism and the probability approach is given by 

\F{t)) = Y J P{n,t)\n) (3) 

m 

It had been shown by Doii that the average of an arbitrary physical quantity B{n) can be 
calculated using the average of the corresponding operator B(t) in according to 

(B(t)) = J2P(n,t)B(n) = (s \ B \ F(t)) (4) 

rii 

with the state function (s \= I- Using the relation (s \ L = the evolution equation for 
an operator B can be written in the form 

d t (B) = (s\[B,L)\F(t)) (5) 
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It should be noted that all dynamical equations govering the classical problem are determined 
by the commutation rules of the underlying operators and the structure of the evolution 
operator L. In our case the dynamics of the model is given by spin-flip processes indicating 
a change of the local spin orientation and through the Kawasaki exchange dynamics. 
The evolution operator for a local single flip-process reads@ 

£/ = £[A(l-d;)4 + 7(l-4H] (6) 

i 

where A and 7 are stateindependent flip rates. A generalization to flip processes with 
restriction is discussed inl!3 

The operators di and d\ fulfil the commutation rule of Pauli-operators. The occupation 
number operator rij = d\di is related to the Ising spin variable by Si — 1 
The evolution operator for spin conserving exchange processes reads 
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Le^EK 1 - + (1 - djfydfijdt] (7) 

1 <i,j> 

with the exchange rate jl. 

The complete dynamics is given by a superposition 

L=pL f + (l-p)L e (8) 

The quantity p represents the probability that the spins follow the Glauber dynamics whereas 
I — p characterizes the amount of spins which are subjected to an exchange process. 
Up to now the flip or exchange processes may be performed independently on the enviroment 
in which the system is embedded and on the static interaction between the spins. The explicit 
coupling to a heat bath and the inclusion of the static interaction can be realized, seet3, by 
replacing the evolution operators Lf, eq.(|^), by 

L f = =u^2[(l- di) exp(-(3H/2)d\ exp((3H/2) 

+ [(1 - 4) exp(-(3H/2)di exp((3H/2)] (9) 

Here, v is a new hopping rate, j3 = T _1 is the inverse temperature of the heat bath and H is 
the Hamiltonian describing the static interaction between the spins. In general the inclusion 
of the mutual interaction leads to nonlocal terms already in the flip evolution operator eq.(^). 
In case of the exchange dynamics the evolution operator L e , eq.(0) should be rewritten in 
the form 

£ e = ^E [(l-^4)exp(-^/2)4^exp(/?i//2) 



+ [(1 - djd\) exp(-(3H/2)d)di exp((3H/2)] (10) 
where /1 is the exchange rate. 
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III. INFINITE RANGE MODEL 



The thermalization introduced by eqs.(|9|,|T(J|), is complete different to the conventional ap- 
proach due to Glauber, seel. The Hamiltonian H, specified below by eg. (|TT|) , mediates the 
coupling to the bath at the fixed temperature T and to the mutual interaction present also in 
the static limit. Physically, the replacement of the operator d\ by exp(— j3H/2)d\ exp(f3H/2) 
in ([D means that a flip-process is realized with a weighting rate exp(/3if/2). After perform- 
ing the flip, manifested by d\, the final state is related to the weighting rate exp(— (3H/2) 
consistent with the fact that only single flip-processes are taken into account. Hence, the 
procedure simulates in an analytical manner Monte Carlo steps. Due to the thermalization 
a spin flip is not independent on the orientation of the other spins. Instead the process is 
self organized by the static coupling of the spins themselves. 

As the simplest case we consider an infinite range model introduced by Kaclll. By reason of 
the long range interaction the model exhibits an exact static solution according to the mean 
field approximation. The Hamiltonian is defined by 

H = -^ E (ii) 

l<i<j<N 

Using this Hamiltonian in the evolution operator eq. (B there appears a conflicting situation 
between the local spin flip process and the infinite range interaction manifested by eq.(|ll|). 
Whereas the long range static interaction favours a reordering of many spins the dynamics 
is refered to single spin flip processes. The situation is comparable with those known from 
spin glasses where the local disorder leads to restrictions within the dynamical processes. 

IV. DYNAMICAL EQUATION 

Using the algebraic properties of Pauli-operators the thermalized evolution operators Lf @ 
is rewritten as: 

L f = z/]T[A(l - di)d]e- E + (1 - d\)d ie E ] 

i 

with E = — — ySj (12) 
2TN 4^ 3 v ; 

3 

Different to the short range model (nearest neighbour coupling) the energy E depends on the 
values of all other spins. It should be noted, that E is obtained if terms with i = j are not 
taken into consideration in eq. (|TTD . When the diagonal terms are regarded the dynamical 
approach does not reproduce the correct stationary state. Within a static consideration the 
mentioned exclusion is irrelevant. 

As the result the flip process is a thermal activated one realized with a certain rate depending 
on the energy E including the temperature. Within the infinite range model it is easy to 
derive that the exchange process is not thermalized. The corresponding terms are canceled 
out mutually. 

Taking into account the analytical form of Lf and L e and using eq.([5|), the evolution equation 
for the order parameter can be obtained: 
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d t (S r ) = up((smhE) - 2(S r cosh E)) + ^(1 - p)J2( S i ~ S r) (13) 

l(r) 

where l(r) means summation over all nearest neighbours of the lattice site r. Note that 
the first part, proportional to p, is originated from the flip process where the second one, 
proportional to 1 — p, can be attributed to the exchange dynamics. The situation is different 
to the conventional analysis where the flip processes contribute to spatial correlations, too. 
The evolution equation for the higher order terms appearing in eq. ([T3|) could be calculated 
in a straightforward manner, seet2. However, in the large iV-limit those terms has not to be 
evaluated provided the system is an ergodic one. It results 

Km E = (E) = ^-(S r ) (14) 

Consequently, eq. flilf ) becomes a closed equation for the order parameter itself. In the 
vicinity of the phase transition the present approach leads to a Ginzburg Landau functional, 
see for instanceil. In a continuous approximation for <p(x,t) = (S r (t)) it reads 

d t ip = vp (sinh(ey?) — 2(p cosh(ey?)) + /i(l — p) VV 
J 2T C 

With 6 = 2T ^ ~T ^ 

The correlation function can be calculated within an expansion in terms of the order pa- 
rameter (S r (t)). In lowest order we find 

CoM)^ . * , ,_ 2 (16) 
The correlation length £ depends on the concentration p and T — T c in the following manner 

The result can be easily interpreted. If the probability p — 1 all spins are subject to 
a stochastic flip-process. But, due to the long range interaction, static fluctuations are 
complete suppressed. The long range force tends to align all spins parallel preventing spatial 
fluctuations. Consequently, the system is not able to establish spatial correlations, hence 
the correlation length is zero. In the opposite case p — > 0, when the whole spins are subject 
to an exchange dynamics, all neighbouring pairs of spins tend to fix their orientation (spin 
conservation). Then, the fluctuations are extremly strong and for p = the system remains 
in a strongly correlated state for finite temperatures. As the result the correlation length 
goes to infinity for each temperature. In the intermediated region < p < 1 the correlation 
length is a product of a p-dependent part and a temperature-dependent part originated from 
the Ginzburg Landau functional. As already emphasized, the situation is complete different 
to the conventional approach where the energy E is a local energy leading to additonal 
spatial fluctuating terms originated from the spin-flip terms. One can also estimate that 
the inclusion of higher order fluctuations does not change the situation. A renormalization 
group procedure based on a Ginzburg Landau expansion starts with the correlation function 
Co obtained in (0). Such an approach yields only a modified critical exponent for the 
correlation length £ m eq . (|17|) . There is also no indication for a tricritical point. 
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V. CONCLUSIONS 



In this paper we have traced back a model with competing Glauber and Kawasaki dynamics 
to a Ginzburg Landau functional where the correlation length depends explicitely on the 
amount of spins which undergo single flip processes and on the temperature. This result 
could be achieved by a coupling to a heat bath and by the consideration of the mutual static 
interaction between the spins. When this interaction is of long range order we obtain a 
closed set of equations for the order parameter. As the result of several conflicting situa- 
tions, long range static force versus a single spin flip process or nearest neighbour exchange 
coupling, or Glauber versus Kawasaki dynamics the correlation length of the system de- 
pends on the probability p which is a measure of the amount of spins which are subject to a 
Glauber process. It should be emphasized that the model is neither model A nor model B in 
the conventional classification^!. Different to this classification the Ginzburg Landau func- 
tional consists of two parts, a nonconserving and spatial independent part and a diffusive 
(conserving) one originated by the exchange dynamics. 
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